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Summary. — We review several aspects of the Sunyaev-Zel'dovich (SZ) effect as- 
sociated with the large scale baryon distribution and its characteristic signatures in 
the statistics of cosmic microwave background (CMB) anisotropies. We discuss (1) 
the contributions to the angular power spectrum from the thermal and kinetic SZ 
effects, (2) the effect of SZ non-Gaussianities on cosmological parameter estimation, 
and (3) the SZ-induced CMB polarization towards galaxy clusters. We also discuss 
the extent to which recent measurements of small angular scale CMB anisotropy can 
be accounted for by SZ clusters and the potential use of SZ polarization towards 
large samples of galaxy clusters as a probe of dark energy parameters. 
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1. Introduction 

The Sunyaev-Zel'dovich effect (SZ; p^) is the upscattering of cosmic microwave back- 
ground (CMB) photons as they pass through the hot gas of a galaxy cluster. This 
upscattering introduces frequency-dependent secondary temperature anisotropics in the 
CMB that are proportional to the integral of the electron pressure H e = n e kBT e along a 
given line of sight through the Universe: 



AT SZ f a T 

(1) = g v {x)y = g v (x) / dryafa) — lie fa) 

J- cmb J m e 



where we have defined y, the Compton y-parameter, <tt is the Thomson cross-section, 
and m c is the electron rest-mass The frequency dependence of the SZ anisotropics 
is accounted for by the function g v {x) — xcoth(|) — 4 with x — ^(A;.bTcmb) > which, 
in the low frequency Rayleigh- Jeans part of the spectrum has the limit g v — > —2. This 
form for the frequency spectrum leads to an apparent temperature decrement, relative 
to the primordial CMB black body spectrum, at frequencies below the null at across, and 
an increment above across- This SZ null occurs at approximately 217 GHz with a slight 
dependence on the electron temperature T e due to higher-order relativistic corrections 
to the inverse-Compton scattering process neglected in the Thomson approximation. In 
Eq. l (Xj l. rj = rj(z) is the conformal distance from the observer at redshift z = to the 
cluster at redshift z, given by 



(2) V (z) 



dz' 
Hiz 7 ) 



where the expansion rate of cold dark matter (CDM) cosmological models with a cosmo- 
logical constant A is 

(3) h 2 = Hi [n m (i + zf + n K (i + zf + n A ] . 

The Hubble distance today is Hq 1 = 2997.9 /i _1 Mpc with h = 0.72 ± 0.08 EJEj. ^ 
denotes the contribution of component i in units of the critical density p cr ;t = 3Hq /8ttG, 
with i = c for the CDM, b for the baryons, and A for the cosmological constant. We 
also define the auxiliary quantities f2 m = fl c + Jlf, and Qk = 1 — fl^ , which represent 
the total matter density and the contribution of spatial curvature to the expansion rate, 
respectively. 



The thermal SZ effect has been imaged in the CMB towards massive galaxy clus- 
ters whose presence is a priori known from optical data 0] E] • The temperature of the 
scattering medium in certain clusters can reach up to 10 keV producing temperature 
changes in the CMB of order 1 mK at Rayleigh-Jeans (RJ) wavelengths. The individual 
galaxy cluster SZ images have a variety of astrophysical and cosmological applications 
including a direct measurement of the angular diameter distance to the cluster through 
a combined analysis with X-ray data and a measurement of the gas mass, and, thus the 
baryon fraction of the universe. We do not discuss these applications of the SZ effect 
in the present review, but rather will concentrate on the statistical study of wide-field 



( ) Note that here and in what follows we adopt units where the speed of light c is set to unity. 
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CMB data where SZ effects lead to anisotropics in the temperature distribution both due 
to resoived and unresolved galaxy clusters. In fact, the thermal SZ contribution is the 
dominant signal beyond the damping tail of the primary anisotropy power spectrum. 

In addition to the thermal SZ effect, the bulk flow of electrons in galaxy clusters and 
other virialized halos that scatter the CMB photons leads to a second contribution to 
temperature fluctuations, the kinetic SZ effect. This effect arises from the well known 
Doppler effect 00: 

cyclop „, (0 

(4) Tf (») = / dr)g(r))h-v g (r),hr)), 

J- CMB Jo 

where v ff is the baryon velocity and 5(77) is the visibility function (also interpreted as the 
probability of scattering within drj of 77) 

(5) g = re- T = X e (z)H T H (l + z) 2 e- T . 

Here t(t)) — J* Q drj' 'f(rf) is the optical depth out to distance ry, X e (z) is the ionization 
fraction as a function of redshift, and 

(6) t h = 0.0691(1 -Y p )Q b h, 

is the optical depth due to Thomson scattering to the Hubble distance today, assuming 
full hydrogen ionization with primordial helium fraction of Y p (= 0.24). From equation 
(£EJ follows the secondary temperature anisotropy due to the kinetic SZ effect. The kinetic 
SZ fluctuation is written as an integral over the product of the radial velocity, n • v g , 
and the density fluctuation associated with the cluster, 5 g : 

^yKinSZ rr]o 

(7) — (n) = / dr)g(r))n-v g (ri : nr))-5 g (r),nri). 

J CMB Jo 

The angular power spectrum of anisotropies generated by the Doppler effect peaks 
at angular scales corresponding to the horizon at the time of scattering. The effect can- 
cels out significantly at scales smaller than the horizon since photons scatter against the 
crests and troughs of the perturbation. The kinetic SZ effect arises from the second-order 
modulation of the Doppler effect by non-linear density fluctuations associated with viri- 
alized objects such as galaxy clusters, and avoids the strong cancellation associated with 
the linear Doppler effect. The kinetic SZ effect is also known as the Ostriker-Vishniac 
effect |S] when the density fluctuations are associated with the linear density field of 
the large-scale structure. In addition to modulation of the velocity field due to den- 
sity perturbations, the velocity field is also modulated by fluctuations in the ionization 
fraction of electrons in a partly reionized universe during the reionization epoch. This 
latter contribution is generally referred to as patchy-reionization [HJ EH ID] • Due to the 
density weighting, the kinetic SZ effect peaks at small angular scales (sub arcminutes for 
ACDM). For a fully ionized universe, contributions are broadly distributed in redshift so 
that the power spectra are moderately dependent on the optical depth r. 

In this review we focus on the statistical characterization of secondary CMB temper- 
ature fluctuations in the form of the angular power spectrum of SZ-induced anisotropies. 
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Fig. 1. - The complex distribution of dark matter (a) found in numerical simulations may be 
replaced with a distribution of dark matter halos (b) with the mass function matching that 
found in simulations. 



SZ maps have important additional applications including number counts of halos where 
a comparison to analytical predictions can be used to derive cosmological parameters. 
These studies will not be the subject of this short review, but are discussed elsewhere in 
these proceedings. 

The outline of this article is as follows. In §2 we derive the power spectrum of 
the thermal SZ effect using the halo model of large scale structure. §3 compares these 
theoretical predictions to the current CMB data. In §4 we derive the power spectrum 
of the kinetic SZ effect. Finally, we discuss secondary CMB polarization due to free 
electrons in galaxy clusters and its cosmological use as a powerful probe of dark energy 
in §5. 

2. Thermal SZ Contribution to the CMB Power Spectrum 

This section reviews the thermal SZ contribution to the angular power spectrum of 
CMB anisotropics following standard derivations published in the literature 1131 1141 

The basic idea of the halo approach to large scale structure is summarized in Fig. 1. 
The complex dark matter distribution is replaced by a population of dark matter halos 
that reproduces the relevant statistical properties of the actual distribution 16 . The 
SZ anisotropy calculation is particularly well suited for the halo approach since high 
resolution simulations of the hydrodynamics of clusters, including radiative cooling and 
feedback from supernovae and galactic winds, indicate that nearly all of the power at 
small scales originates from regions associated with virialized halos rather than from fila- 
mentary structures or the diffuse gas between halos JS]- The semi- analytic calculations 
of Rcfs. ^3 E| are extended here by accounting for the self-gravity of the cluster gas 
and allowing for scatter in the concentration-mass distribution. 

To calculate the angular power spectrum associated with the thermal SZ effect, it is 
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convenient to describe statistics in terms of the baryon gas pressure Tl g , which is related 
to the pressure of the electrons in the cluster by il e (r) = |^±|n s (r), where X = 0.76 is 
the primordial hydrogen abundance. We consider the difference in the local gas pressure 
relative to the mean pressure averaged over all positions, at a given epoch, and write the 
fluctuation as <5n(x, t) = (II — ft)/n, where II denotes II e or n g . The mean value of Su 
vanishes by definition: 

(8) (Sn) = Hm J|x|< f. p3 = . 

The correlation function of pressure fluctuations is defined as £(x) = (<5n(xi)<5n(x 2 )), 
where homogeneity and isotropy imply that correlations only depend on the absolute 
separation of the two locations, x = |x x — x 2 |. If the fluctuations are Gaussian, i.e. the 
random field related to Sji has a Gaussian probability distribution, then the distribution 
function is fully specified by the correlation function alone. A generic prediction of in- 
flationary cosmology is the Gaussianity of the spectrum of initial fluctuations. However, 
non-linear growth and the formation of virialized halos lead to significant non-Gaussianity 
in the pressure field at the low redshifts probed by the SZ effect. Hence, to fully spec- 
ify the statistics of the SZ power spectrum one must make measurements beyond the 
correlation function. 

Instead of the real space correlation function, we will discuss statistics primarily in 
Fourier space where the fluctuation in pressure is 



(9) <5 n (k) = / d 3 x<5 n (x)e 



-ik-x 



It is easily shown that Gaussianity of <5n(x) implies Gaussianity of 5n(k). We define the 
three-dimensional power spectrum of pressure fluctuations in the large scale structure as 

(10) (Mk)Mk')) = (2nfS D (k - k')Pnn(fc) . 

The Dirac delta function 5d expresses the fact that wave modes k and k' are independent 
and uncorrelated since the different (5n(k) have random phases. Spatial isotropy requires 
the power spectrum to be independent of the phase of the wave vector and depend only 
on the magnitude of the vector. The power spectrum and the correlation function are 
related via 



(11) Pnn(fc) = J d 3 x(,(x)e-^ = J dx4nx 2 £_(x) 



sin(fca;) 
kx 



The halo model for large scale structure predicts two separate contributions to the 
power spectrum: 

(12) Pnn(k;z) = P^ 1 (k;z)+P^ l (k;z), 
where 

P^{k-z)=ll nu {k-z), and 

(13) P^{k-z) = [ll nn {k-z)fP{k;z), 
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Fig. 2. - Left: The density (top) and pressure (middle) profiles for the self-consistent self- 
gravitational model (solid) and the isothermal cluster model (dotted) compared with numerical 
simulations from the Santa Barbara cluster comparison project. The isothermal model over- 
estimates the pressure and density in the inner and outer regions of the cluster. The lower 
panel shows the ratio of the pressure profile calculated with self-gravity to the profile calculated 
without self-gravity for /b = 1/10 (solid), the value used in the Santa Barbara cluster compari- 
son project, and fb = 1/6 (long-dashed), our fiducial value. Self-gravity raises the pressure by 
~ 20% outside the core for our fiducial value. Right: The mass dependence on the dark matter 
power spectrum (a) and the pressure power spectrum (b). Here, we show the total contribution 
broken up into mass limits as labeled in the figure. As shown in (a), the large scale contribution 
to the dark matter power comes from massive halos while small mass halos contribute at small 
scales. For the pressure, in (b), only massive halos above a mass of 10 14 Mq contribute to the 
power. 



denote the power spectrum for correlations of gas pressure within a single halo (lh) 
and between two separate halos (2h). It is found that the 1-halo term dominates on 
small angular scales. The functions 1° n n an d ^i".nn are defined below and P(k; z) is the 
linear matter power spectrum as a function of redshift z. Useful fitting formulae for the 
transfer function relating P(k; z) to the primordial power spectrum P{k) may be found 
in Ref. 

Using the three-dimensional Fourier transform of the (spherically symmetric) halo 
pressure profile II(r), 

(14) n e (fc) = j d 3 rn e (r)e- ikx = J°° dr 4nr 2 n e (r) S -^j^ , 
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we can write the two /-terms in Eq. I|13|) as 

7 2 ° nn (fc; z) = J dM J dc -^j c |n e (fc, M, c; z)\ 2 and 
f f d 2 n 

(15) Il nu (k;z) = J dM J d c ^-6 1 (M;z)n e (fc,M,c;z). 

Here, d 2 n/dMdc is the bivariate halo mass function, the co-moving number density n 
of halos between mass M and M + dM with concentration parameter between c and 
c + dc. The concentration parameter c = r v /r s is defined as the ratio of the virial 
radius of the cluster r v and a characteristic scale radius r s . In writing Eq. we 
assumed that halos trace the linear density field such that the halo power spectrum is 
Phh(k; z\M) = b\(M)P(k] z) where M denotes the halo mass and the halo bias b\(M\ z) 
is 113 ' 

(16) bl (M; Z ) = l + a ^ M ^- 1 ■ 



S c (z) 6 c (z){l + [av 2 (M;z)]p}' 

where (a,p) are numerical parameters and u(M; z) = S c (z)/a(M; z) as discussed in more 
detail in section 2.1. We see from Eq. I|15|) that the two ingredients involved in calcu- 
lating the power spectrum with the halo model are the radial pressure profile of the gas 
within each halo and the mass function, which specifies the number of such halos. 

The angular power spectrum of the thermal SZ effect is just proportional to a projec- 
tion of the three dimensional pressure power spectrum integrated along the line of sight. 
To derive the angular spectrum, we take spherical harmonics of Eq. (1) such that 

(17) ^W-Ef^) Y ^)- 



irn 



Making use of the Rayleigh expansion of a plane wave 

(18) e* k -" = 4nY,*\n(kvW rn (m m (n) , 

Irn 

we find the spherical harmonic moments of the SZ map 

(19) (^T) =9u^ w| ^n e (k)Y4(k) I dna-^Mkr,). 

The angular power spectrum of the SZ map, is defined in terms of multipole moments 
[AT/TV as 

(20) ( [AT/ T] i m [AT/ T] *, m , ) = Cf z 6 w S rnm , . 

Using equation (0 and the relation (n e (k)n e (k')) = fl^ [(<5 n (k)<5 n (k')) + 1] we obtain 



Cf z = g l(^L) [n e {Q)f- fk'dkPun 
\m e J 7T J 



(fc)x 
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(21) X JdViJ dr, 2 a- 2 {i 1 i)a-' 2 {T l2 )ji{kr ll )ji{kr]2), 

where n e (0) is the mean electron density today and we used fi e (r]) = n e (0)a(r])~ 3 . 

With the Limber approximation |21| . or, equivalently using an approximate identity 
related to the completeness of spherical Bessel functions 



(22) / dfcfc 2 F(fc)j ; (fc77)ji(W) « ^ 2 S D ( V - V ')F(k)\ k=l/dA , 
the projected SZ power spectrum can be written as 

(23) Cf z = fdr ^j^Pnu ( k = -j~rT '' z ( r ) 



d\{r) V" d A {rY 

where d A is the comoving angular-diameter distance, which in the case of a flat-cosmological 
model is d A = rj. We have defined the SZ weight function 

The Limber approximation assumes that F(k) is a slowly varying function when com- 
pared to oscillations in ji (x). Under this assumption, contributions to the power spectrum 
arise only from correlations on equal time slices of spacctimc. 

Given that the three-dimensional power spectrum of pressure can be written as a 
combination of 1- and 2-halo terms and the fact that the 1-halo term dominates the 
fluctuations, we can simplify the integral further and write it as |15| 

where d 2 V/dzdVL is the physical volume per unit redshift per unit solid angle and 11/ is 
related to the three-dimensional Fourier transform of the electron pressure profile via 

w sz (z)n e (k = -J-. 

(26) H(z) = i - 



Using the projected angular size of a cluster 6 = r/d A , one can simplify further and 
reduce LI/ (z) to the following two-dimensional integral 

/■flout r / 1 \ 1 f+'out 

(27) I1 1 {z) = 2itW sz {z) j WJ [(l + -)o\J U.(r\\,r ± = d A 0)dr\\ , 

where we have rewritten the volume integral Eq. 1141 in terms of two separate integrals 
involving the line of sight distance through an individual halo, rii , and the 2-dimcnsional 
projected extent of the cluster, rj_ = d A 8- Here, r out is the outer radius of the halo, 
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Fig. 3. - The angular power spectra associated with the thermal and kinetic SZ effects. As shown 
in (a) , the thermal SZ effect is dominated by individual halos, and thus, by the single halo term 
(lh, dotted line), while the kinetic effect is dominated by the large scale structure correlations 
depicted by the 2- halo term (2h, dashed line). In (b), we show the mass dependence of the 
thermal and kinetic SZ effects with a maximum mass of 10 16 and 10 13 M . The thermal SZ 
effect (solid lines) is strongly dependent on the maximum mass, while the kinetic effect (dotted 
lines) is not since it receives most of its contribution from large scale correlations. 



which with a pressure profile that falls off with distance can be taken as r out — > oo, while 
#out = r OVLt /d,A- 

We now describe the ingredients required for the halo-based calculation of the SZ 
power spectrum in detail. 

2T. Halo Mass Function. - The mass function of the halo distribution is derived 
from the Press-Schechter formalism [22]. Press-Schechter theory is based on the idea 
that fluctuations in the linear density field with S > 5 C detach from the local Hubble 
expansion of the universe and collapse to form non-linear structure. The prediction for 
the fraction of the volume that has collapsed is 

9 /*°° 

(28) f coli (M(R),z)= / dSe-t'WM, 

V 2iT(j(R, z) Js c 

where R is the radius over which the density field has been smoothed, which is related to 
the halo mass by M(R) = p m ^rR 3 with p m the co- moving matter density of the universe. 
The number density of halos is then found to be given by [221 12H] 

(2q s dn(M, z) _ p m df coU (M(R), z) _ p m dv 

[ ' dM ~ M dM ~ M'^'dM' 
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where 



2A 2 



2 



(30) /(l/ ) = A /r^L[i + ( OI/ 2)-P] e - 
Here 

(31) i/(M;*) = -^L 

cr(M;z) 



and 



r / x 3(12tt)3 
(32) S c (z) ~ -^i- 



5 , / ! <>. 



1 - _ln 1 



936 V n m {l + z) 3 



is the critical density required for spherical collapse at a redshift z in a flat ACDM model 
(see e.g. Ref. [21] and references therein for this and related formulae). The variance in 
the density field smoothed with a top-hat filter of radius R = (ZM / An p m )$ is 

(33) <J 2 (M; z) = G 2 (z) / ^^V(fci?)| 2 , 



k 2tt 2 



where 



3 

(34) W(x) = -yfsinfx) - xcos(x)] , 

x 6 

P{k) is the linear matter power spectrum, and 

H(z)f™dz'(l + z>)[H(z')}- 3 



(35) G{z) 



H ^dz'{l + z')[H{z')]- 3 



is the linear theory growth factor (often also denoted D(z)). A, p, and a in equation l|3UI) 
are constants, with the canonical Press-Schechter (PS) and Sheth-Tormen (ST) mass 
functions corresponding to the parameters (p — 0,a — 1) and (p — 0.3, a = 0.707), 
respectively. The normalization A is determined by requiring mass conservation such 
that 

(36) — / dMM^- = / d„/(„) = 1 . 



p m J dM 

2'2. Dark Matter Density Profile. - The halo mass function has to be supplemented 
by the dark matter density profile. We define the dimensionless variable u = r/r s , where 
r s is a characteristic scale radius. The dark matter profile within each halo is then defined 
as an NFW profile gS] 

(37) p d (r) = p a Qd(u) = p s w _1 (l + u)~ 2 . 

Within the context of the spherical collapse model, the outer extent of the cluster is 
taken to be the virial radius 



(38) 



3M 



47rp m (z)A(z) 
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where p m (z) oc (1 + z) 3 is the mean background matter density of the universe at redshift 
Z, and 



(39) 



A(z) ~ 18tt 2 



215 V^m(l + ^) 3 



is the overdensity of the halo relative to the background density [21] • The ratio of the 
virial radius to the scale radius is called the concentration parameter c = r v /r s . Together, 
c and M completely determine the dark matter distribution of a given halo. 

2'3. Concentration-Mass Distribution. - Numerical simulations indicate that clusters 
of a given halo mass have a range of concentration parameters. To account for this 
distribution we describe the number density of halos of mass M and concentration c by 
the bivariate distribution function, 



(40) 



dMdc 



(M, c; z) 



dn 
dM 



(M;z)V(c\M;z). 



where V{c\M\ z) is the probability of a halo having concentration c given that it has mass 
M. Numerical simulations 26 find an approximately log-normal distribution 



(41) 



T{c\M-z)dc 



exp 



[log C -logg(Af;z)] ; 



dc 



2na\ 



ln(10)c ' 



where the mean concentration parameter c is related to the halo mass via 



(42) 



c(M ; z) 



co 



M 



M*(z = 0) 



M*(z) is the mass scale at which i/(M*;z) = 1 and Cq and a c are constants whose 
numerical values are typically chosen to be cq — 9 and a c = 0.13. It is in general possible 
that ciogc = c\ogc(M', z ), a s suggested by the scaling arguments of Ref. [27j and indeed a 
mass dependence has been noted in numerical simulations |26| . However, for simplicity 
we fix the width of the concentration distribution to be <7i og c = 0.14 independent of mass. 

2'4. Gas Pressure Profile. - The gas density profile in terms of the dimensionless 
parameter u is pg(r) = poQ g (u). The isothermal assumption for the gas distribution 
used in Refs. |13l I14j was found to be inconsistent with the properties of numerically 
simulated clusters |17| . We therefore implement a polytropic equation of state and write 
the gas pressure as 



(43) 



P g (r) = P oP (u) = P [q s (u)]' 



with g g (0) = p(0) = 1, so that po and Pq correspond to the central density and pressure. 
The gas profile within each halo has to self-consistcntly satisfy hydrostatic equilibrium, 



(44) 



1 dP n 



GM(r) 



Pa 



dr 
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where M(r) = dr' '47rr /2 [p g (r') + pd(r')] is the total mass of baryonic gas and dark 
matter. Rewriting Eq. (|44|l in terms of the dimensionless variables we introduced above 
leads to the following second order differential equation 



(45) 



d_ 

du 



2 T 



du 



= -\u 2 [g d (u) + /3g g (u)} 



where 

(46) 

and 



(3 



Po 

Ps 



(47) 



A 



Po 



In general g g = g g (u; A, /3, 7, c, /b). To fix g g (u; c, /b) for a given halo of concentration 
c and cosmological baryon fraction Jb = Slb/^m we impose the boundary condition 
that the gas density profile traces the dark matter density profile in the outer region of 
the halo. This ansatz from Ref. \Tj\ ensures that the hydrodynamical properties of the 
cluster resemble those found in simulations of clusters which include both gas and dark 
matter. In order to implement this boundary condition and determine A(c; /b), /3(c; /b), 
and 7(0; /b), we minimize the functional 



(48) 



V[g s (u;\,j3,j,c,f b )] 



— hsM ~ s d (u)} 2 
1 „ u 



subject to the normalization constraint 

p g (r v ) f2 b 



(49) 



Pd(^v) 



/3c(l + c) 2 e g (c), 



where s(u) = dlng/dhiu is the local logarithmic slope of the density profile g(u). We find 
that this procedure matches that of Ref. when the gas self-gravity is neglected. 

In Fig. [21 (left panel) , we present a brief comparison of our model and the isothermal 
model with results from numerical simulations done as part of the Santa Barbara cluster 
comparison project [2HI - As shown in Fig. |2 the isothermal cluster overestimates the 
gas density and pressure in the inner region and also shows a departure at the outer 
extent of the cluster. Comparing our model with that of Ref. ^Kj we see that including 
self-gravitational effects raises the pressure by ~ 20% outside the cluster core, which can 
increase the amplitude of the angular power spectrum of the SZ effect by as much as 
25-30%. 



3. Comparison with Current CMB Data 

In 2002, two significant results on CMB temperature fluctuations at arcminute an- 
gular scales appeared. The Cosmic Background Imager (CBI; j2Hl) and the Berkeley- 
Illinois-Maryland Array (BIMA; POj) measured an excess of small-scale power beyond 
the damping tail of the primary CMB. Interpreting the CBI and BIMA results as low 



STATISTICAL IMPRINTS OF SZ EFFECTS IN THE COSMIC MICROWAVE BACKGROUND 



13 




Fig. 4. - (a) Small-scale CMB anisotropies as observed by the CBI and BIMA experiments. The 
three solid lines show predictions for the SZ effect for varying erg: 1.0 (upper), 0.9 (middle), 
and 0.8 (lower). For reference we show the predicted CMB primary anisotropies for our fiducial 
model. The solid error boxes show the published Gaussian errors, while the dotted boxes 
include the non-Gaussian trispectrum. When comparing data to predictions, we account for 
the full window function of these observations, (b) Constraints on cosmological parameters. 
We consider constraints on n m and as using small-scale power and show 68%, 95%, and 99% 
confidence limits. 



redshift SZ contributions to the temperature anisotropies allows constraints to be placed 
on cosmological parameters. The cosmological implications of the small-scale excess in 
CBI were first studied in Ref. [2] using numerical simulations of the SZ effect associated 
with unresolved clusters. A more detailed analysis of the small-scale CBI data, based on 
semi-analytic calculations, was presented in Ref. |15) . where a constraint was placed on 
the normalization of the matter power spectrum erg. 

To study the constraints on cosmological parameters from small-scale anisotropies, one 
generally performs a likelihood analysis. Here, one must take into account the window 
functions of band-power from CBI and BIMA. In order to compare data to theoretical 
predictions, we calculate the relative x 2 between data and models 

(so) x 2 = £(Ci-ei)<V(^- c i)> 

i<j 

where Ci is a band-power estimate from recent observations, Ci is the model prediction, 
and is the inverse covariance matrix. As we will now discuss, due to the highly 
non-linear nature of galaxy clustering, one must take non-Gaussian contributions to the 
covariance into account. 

3'1. Covariance. - The measurements of the power spectrum reported by the CBI 
and BIMA experiments are binned estimates of power in multipole space with window 
functions W(lj) at each band i, where 1^ is a vector in a plane of the sky sufficiently small 
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to be considered fiat. One can write these band-power estimates as 



(51) 



c ; = I ^w(h)w(-h)e(i t )e(-h) , 



where Ai — J d 2 \iW 2 (\i) is the window-weighted area of the two-dimensional shell in 
Fourier space and 0(li) is the Fourier transform of temperature fluctuations O(n) = 
AT/T(n) in the SZ map. Following Ref. [22, we can write the covariance matrix of 
these band-power measurements as 



(52) 



Cij — 



(2tt) 

A 



A, + (2tt) 2 



where 



(53) 



rpS Z 



A,, 



and A is the total area of the survey in steradians. The first term in equation l|52l) 
is the usual Gaussian sample variance and the second term is the extra non-Gaussian 
contribution that arises from the trispectrum associated with the SZ effect. The Gaussian 
noise also includes an additional noise variance due to instrumental noise. In the case 
of CBI and BIMA, the first term, including instrumental noise, is generally determined 
during the band-power measurements. Thus, we only consider the extra covariance from 
the non-Gaussian contribution. We again use the halo approach to model this term. At 
the small scales relevant for current data, the one-halo term dominates We thus 

consider the trispectrum due to single halos, 



(54) T sz (h,\ 3 )=T sz (h,-h,\^-l 3 



( h 



dr W^rf 

d\{r) UU \d A (ry d A (r) 



;z(r) 



Note that the projected trispectrum is generally a function of four 1 vectors and may 
therefore be represented by a general quadrilateral in Fourier space. In the case of the 
power spectrum covariance we are only interested in terms involving Tgz(lj, — lj, lj, 
i.e. parallelograms which are defined by either the length Uj = |li — lj| or the angle 
between 1^ and Ij. This can be noted from the fact that the correlator between (CiCj) 
involves the temperature fluctuations as (9(li)9(— li)9(lj)9(— lj)) which is simply equal 
to the trispectrum with the above dependence. 

Under the halo model the trispectrum involves four terms related to configurations 
with 1 to 4 halos. On small scales the 1-halo term dominates, similar to the case of the 
power spectrum, and we can write this contribution to the trispectrum responsible for 
power spectrum covariance as 



(55) 



7^(k 1 ,k 2 ;2) = J dM J dc-^- c |H(|ki |, M,c;z 



)| 2 |n(|k 2 |, m,c A z)Y 



In addition to the extra variance (when |li| = |lj|), this non-Gaussian term also con- 
tributes to the covariance and correlates band-power measurements between different 
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bins (when |lj| ^ |L,|). Note that in the case of the 1-halo term, the trispectrum is 
independent of the angle between 1^ and lj and only depends on the amplitude of these 
two vectors, though, this is not the case in general when other terms, such as 2-halo and 
higher, are involved. 

In Fig. 4, we show the extra variance contribution to the band-power estimates due 
to non-Gaussianitics. These non-Gaussian errors are such that they decrease with in- 
creasing sky area of the survey and should be understood as arising from the rareness 
of galaxy clusters that contribute to the SZ anisotropics. In Fig.0] we show the angular 
power spectrum of the SZ effect along with the primordial anisotropics for our fiducial 
cosmological model. Here we plot flat-band power per logarithmic interval in /xK such 
that 



(56) AT= ^1+J1c^Tcmb, 

where Icmb = 2.726 x 10 6 /iK. We also show flat-band power estimates of the anisotropics 
at small angular scales from CBI J55| and BIMA j3U] ■ The two error boxes represent the 
Gaussian error published by both experimental groups and the total which includes 
the non-Gaussian error. For comparison, we show three predictions for the SZ power 
spectrum for as ranging from 0.9 to 1.1. In general the SZ angular power spectrum Cf z 
scales approximately as (as) 7 

To extract cosmology, we consider a parameter set comprised only of (f2 m , as), fixing 
all other cosmological parameters to fiducial values given by the WMAP best-fit |3j • The 
1-a and 2-a constraints on this two-parameter space are shown in Fig.0] Our constraints 
are consistent with those obtained by Ref. JS] and suggested in Ref. [3T] . Marginalizing 
over values of Q m in the range of 0.14 to 0.96, we find a 8 {Q h h/0.036) - 3 = 0.95±°;2° at 
95% confidence. The data allow as values of around 1.0 to 0.85 at the la level, consistent 
with complementary cosmological measurements. The inclusion of self-gravity and the 
distribution in concentration increases power beyond the models in Ref. ^7]. Note that 
our results, while consistent with results from the analysis in Ref. |84| . still allow a larger 
range of values for as than suggested there. 

While small scale anisotropics can be interpreted as being due to the thermal SZ effect, 
this interpretation is not unique. Arcminute-scale anisotropies can easily be generated 
at the surface of last scattering if the primordial power spectrum had a feature. It 
is, however, important to keep in mind that scale-dependent features in the primordial 
power spectrum may be unnatural in specific models of inflation. Other effects producing 
small-scale anisotropies include primordial voids at last scattering. If the rcionization 
optical depth is high, as suggested by the recent WMAP results and the reionization 
process is associated with massive stars, then a substantial small-scale signal could also 
be generated by the thermal SZ effect associated with the first supernova bubbles |35j . 
Similarly, at low redshifts, the secondary signal could be due to foreground sources such 
as unresolved point sources. 

As discussed in Ref. |36| , to uniquely determine if the contribution is associated with 
the thermal SZ effect, one needs to establish the frequency spectrum with multi-frequency 
observations, or consider cross-correlation studies with tracers of large scale structure. 
If the contribution is generated at last scattering, then one does not expect a significant 
correlation with, say, the low redshift galaxy distribution, while a strong correlation is 
expected, if arcminute scale fluctuations are related to the thermal SZ effect. A similar 
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cross-correlation analysis would be required to separate the thermal SZ effect in galaxy 
clusters from the effect related to the first supernovae. The latter signal would not be 
resolved, while with high resolution SZ maps, the removal of resolved SZ clusters will 
lead to a substantial decrease in the amplitude of the SZ power spectrum, by which the 
presence of SZ supernovae can be explored. 

In addition to the power spectrum, due to the highly non-linear nature of the thermal 
SZ effect, a significant non-Gaussianity is induced by SZ contributions to arcminutc- 
scale CMB fluctuations El IH which can be measured from the bispectrum, the 
three-point correlation function in Fourier space, or more easily, using measurements 
of the skewness, among others. The SZ effect is also strongly correlated with angular 
deflections associated with gravitational lensing modifications of CMB anisotropics 
1371 |2U • Future experiments such as Planck will be able to measure some of these higher- 
order effects. 



4. Kinetic SZ Contribution to the CMB Power Spectrum 

We next provide a brief discussion of the kinetic SZ contribution to the CMB power 
spectrum. The kinetic SZ effect results from the modulation of the large-scale velocity 
field by non-linear density structures such as galaxy clusters. In Fourier space, the effect 
is described as a convolution of the density and the velocity fluctuations (cf. Eq. J7J). If 
we consider the flat-sky limit in which the Limber approximation (21) applies, then we 
can write the associated power spectrum as 

Recall that 17(77) and G{rj) are the visibility function © and the linear growth function 
((35(1 . respectively. The mode-coupling integral is given by 

(58) I v (k) = / dk t f +1 dn {1 '^ )il 2 ' 2fiyi) P88(ky 1 )P gg (ky 2 ) , 

J J-i 2/2 

where P55 and P gg denote the power spectra of perturbations in the dark matter and the 
cluster gas densities, respectively. In the above, fj, = k • k\, y\ = k\jk and 7/2 = k 2 /k = 
yj\ — 2fiyi + y\. The wave vectors ki and k2 = k — ki capture the convolution between 
velocity and gas density perturbations. We refer the reader to El CUB EH E] f° r the 
details of the derivation of equations 1(57)1 and ifSS)! . The velocity field power spectrum is 
related to the linear dark matter density field through the use of linear theory arguments 
involving the continuity equation where v = — iGSh/k 2 , which was used in deriving 
Eq. P| . 

Since velocity fluctuations have a much larger coherence length than the non-linear 
gas density fluctuations traced by galaxy clusters, we consider the limit in which the 
velocity field is uncorrelated with the clusters. Physically, this can thought of as the 
limit in which large-scale bulk flows are modulated by small-scale overdensities. In this 
small-scale limit, the dominant contributions to the integrals in equations (|57(l and 1(58(1 
arise from modes with k 2 — |k — ki| ~ k such that yx <C 1 and y% — > 1, and equation 
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Fig. 5. - Line of sight projected maps of the thermal (left) and kinetic (right) SZ effects (Figures 
reproduced from Ref. }42j ) . The maps are 1° x 1° and cover the same field of view. Note that 
the thermal SZ map picks out massive halos, while contributions to the kinetic SZ effect come 
from a wide range of masses. Unlike the thermal SZ effect, which produces a negative decrement 
at Ray leigh- Jeans wavelengths, the kinetic SZ effect varies between negative and positive values 
depending on the direction of the velocity field. Here, structures in red are moving towards the 
observer while those in blue are moving away. 



(|57|) reduces to 

(59) Cf* SZ = i / d V { -^P gg (k)v^ s , 
where v^ ma is the rms of the (large scale) uniform bulk velocity 

(60) , r 2 ms = J dk^l . 

The factor of 1/3 in Eq. arises from the fact that the rms in each component is 1/3 
of the total velocity. 

Just as in the halo model calculation of the thermal SZ effect, we split the power 
spectrum of perturbations in the cluster gas densities into 1-halo and 2-halo terms, P gg = 
Pgg + Pgg - For the thermal SZ effect we found that the 1-halo term was the dominant 
contribution and that correlations between halos could be neglected. For the kinetic 
SZ effect, however, we find that the 2-halo correlation term is most important and the 
1-halo term can be neglected. This is illustrated in Fig. [31 where we show our prediction 
for the CMB power spectrum associated with the kinetic SZ effect and a comparison 
with the thermal SZ contribution. It is observed that the kinetic SZ contribution is 
roughly an order of magnitude smaller than the thermal SZ contribution. We also see 
that the thermal SZ effect is dominated by individual halos, while the kinetic SZ effect 
is dominated by large scale structure correlations depicted by the 2-halo term. We 
understand these results as follows: The temperature anisotropies associated with the 
thermal SZ effect are proportional to the temperature weighted electron density (or the 
electron pressure, cf. equation (1)) and are therefore dominated by the most massive 
clusters with a correspondingly high electron temperature. The kinetic SZ contribution 
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is independent of the cluster temperature and depends only on the electron density and 
the cluster velocity. Contributions to the kinetic SZ effect therefore come from baryons 
tracing all scales down to small mass halos and the CMB power spectrum associated 
with the kinetic SZ rescattering is dominated by the large scale correlations of the halos. 

The different mass dependence of the two effects suggests that wide-field thermal and 
kinetic SZ maps will show characteristic differences in that massive halos, or clusters, will 
be clearly visible in a map tracing the thermal SZ effect, while the large scale structure 
will be more evident in a kinetic SZ map. As shown in the thermal and kinetic SZ maps 
in Fig.[S]from |42| . numerical simulations are in fact consistent with this picture (see also 

ESI)- 

In addition to the contribution due to the line of sight motion of halos, there is an 
effect resulting from halo rotations as discussed by 03]. Recent studies have considered 
this non-uniform motion of cluster gas with numerical simulations (e.g. |45p. The signal 
is considerably smaller than the bulk flow kinetic SZ effect and therefore unlikely to 
be of interest to current and next-generation measurements of arcminute-scale CMB 
anisotropies. 

5. SZ Contributions to the Polarized CMB 

Polarization of CMB anisotropies is generated when the CMB photons scatter off free 
electrons. The polarization therefore traces the ionization history of the universe. The 
universe was fully ionized at early times, before last scattering occurred at a redshift of 
about z rec f 100, after which the universe became neutral and radiation and matter 
decoupled. Last scattering generated the primary polarization of the microwave back- 
ground radiation. Since the generation of polarization is a strictly causal process this 
signal is expected to peak at the horizon scale at recombination. Causality doesn't allow 
a polarization signal on larger scales. Such a signal, however, has recently been measured 
by the Wilkinson Microwave Anisotropy Probe (WMAP) ^| . This large-scale secondary 
polarization is interpreted as a signature of a late time reionization of the universe at a 
redshift of z ~ 10—20. Since reionization leads to free electrons in galaxy clusters we also 
expect small-scale secondary polarization |47l I48| . In Ref. |49| . we revisited the problem 
of measuring a CMB-induced polarization signal towards resolved galaxy clusters. This 
was previously studied by Refs. [Tl 1511 1521 I53| . 

Linear polarization of the cosmic microwave background is generated through rescat- 
tering of the temperature quadrupole. In a cosmological model with dark energy the 
quadrupole evolves between the last scattering surface (z — 1100) and us (z = 0) due the 
integrated Sachs- Wolfe (ISW) effect. The quadrupole- induced polarization signal there- 
fore provides an opportunity to probe dark energy through the ISW effect. Kamionkowski 
and Loeb |541 15^) have considered the possibility of using multiple such measurements 
to reduce the cosmic variance uncertainty in the CMB temperature quadrupole. The 
connection between properties of dark energy, the quadrupole at the cluster redshift and 
polarization in the direction of galaxy clusters is summarized here. 

5T. CMB-induced Polarization towards Clusters. - CMB polarization towards clus- 
ters is generated when the incident radiation has a non-zero temperature quadrupole 
moment. The two dominant origins for this quadrupole moment are: (a) a projection 
of the primordial CMB quadrupole to the cluster location, and (b) a local kinematic 
quadrupole from cluster peculiar motion. Towards a sufficiently large sample of galaxy 
clusters, we can write the total rms degree of polarization as f^ota! = Pp Ilm + Pxim 
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Fig. 6. - The temperature quadrupole, Cf e , (O = AT/T), as a function of redshift. We show 
both the primordial and the kinematic quadrupole. The bottom kinematic-quadrupole curve is 
for g(x) = 1, as appropriate for the Rayleigh- Jeans (RJ) part of the frequency spectrum, and 
the dashed and dotted curves are, respectively, for frequencies 150 and 220 GHz. 



Ref. [51] . where 

(61) Fp r im oc (r) Q rms (z) (X 7C^) , 



(62) P KiD cx g(x) (t) (ft 2 ) . 

T = a T J d£n e (£) is the scattering optical depth of a cluster with I the physical line 
of sight distance through the cluster and n e (£) the electron number density. Since 
we are averaging over large samples of clusters, we consider the sample-averaged op- 
tical depth, (r). (3 t = v t /c gives the transverse component of the cluster velocity and 
g(x) — {x/2) coth(x/2), with x = hv/ksTcMB, is the frequency dependence of the kine- 
matic effect. With the optical depth in individual clusters determined by complementary 
methods, such as the Sunyaev-Zel'dovich (SZ, |T]) effects, one can invert the measured 
polarization, equation l(HT|l . to obtain the rms CMB-quadrupole, Q rms (z) = (5C2/47T) 1 / 2 , 
at the cluster redshift. 

5'2. Primordial Quadrupole. - The primordial CMB-quadrupole is dominated by two 
effects. The Sachs- Wolfe (SW, jSS]) effect arises as a combination of gravitational redshift 
and time-dilation effects and can be viewed as a direct projection of the conditions at 
last scattering with no evolution after that time: 

AT\ SW _ *fa = »fc) 
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Fig. 7. - CMB polarization due to galaxy clusters. The polarization vectors give a representation 
of the expected polarization from a cluster at the corresponding location in the sky. The scale 
is such that the maximum length of a line corresponds to a polarization of 4.9r /iK, where r is 
the optical depth of the cluster. The background color represents the temperature quadrupole. 
The left map shows the resulting polarization contribution due to scattering of the primordial 
temperature quadrupole alone. The map for the total polarization contribution (primordial and 
kinematic) is indistinguishable from this map. To make the kinematic polarization visible, we 
arbitrarily increase its amplitude by a factor of 100 (right map). The primordial polarization 
dominates the total contribution even at high frequencies where the kinematic quadrupole is 
boosted due to its spectral dependence. 



where $ is the Newtonian potential. The integrated Sachs- Wolfe (ISW) contribution 
arises along the photon path from the time of last scattering to today, as the CMB 
photons pass through a time-varying potential: 

/AT\ ISW f m 
(64) (-) =2|_«„. 

Effectively, the photon receives a shift in energy because the potential it falls into is 
different from the potential it must climb out of. The ISW effect is absent in a matter- 
dominated, critical-density universe (Einstein-de Sitter). In a universe with dark energy 
(w = p/p < 0) or a cosmological constant, A, (w = —1) the ISW effect leads to an 
increase in power on large scales. The expected redshift evolution of the quadrupole, 
Ci=2(z) — Cf^(z) + Cj^(z), is hence characterized by a rise at low redshifts (z < 1), 
the time at which the universe becomes dark energy dominated. 

We calculate C 2 ee (z), 9 = AT/T, following Ref. g7] for our fiducial ACDM cosmol- 
ogy and show the result in Fig. El Note that at redshift z = 0, with the power-spectrum 
tilt fixed at unity, COBE finds C 2 ee (z = 0) = (27.5 ± 2.4/iK) 2 [Hlj. At higher redshifts, 
the mean quadrupole moment evolves due to the integrated Sachs- Wolfe effect and pos- 
sibly, if the power spectrum is not flat, due to any scale dependence since the quadrupole 
probes smaller distances at earlier times. The primordial quadrupole has a coherence 
length comparable to the horizon. We thus expect the polarization signal measured on 
a O(10°) patch of sky may differ by order unity from that on a different 0(10°) patch 
of sky. Finally, note that Thomson scattering from cold electrons will not change the 
photon frequency. Thus, there will be no frequency dependence of the rescattering of the 
primordial quadrupole. 
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5'3. Kinematic Quadrupole. - The origin of the kinematic polarization effect is under- 
stood as follows. Consider electrons moving with peculiar velocity, fi = v/c, relative to 
the rest frame defined by the CMB. The Doppler-shifted spectral intensity of the CMB 
in the mean electron rest frame is 



(65) 



C 



e xj{l+0iJ.) _ i ' 



where x = hv/ksTcMB, 7 = (1 — /3 2 ) -1 / 2 and \i is the cosine of the angle between the 
cluster velocity and the direction of the incident CMB photon. When expanded in terms 
of Legendre polynomials, the intensity distribution is 



(66) 



C- 



e x - 1 



I Q + Ix/J, 



2(e x -\) 2 



x z fi z [l 



which contains the necessary quadrupole component under which scattering generates 
polarization. 

Unlike the primordial quadrupole, the kinematic quadrupole has a frequency depen- 
dence which we denote by g(x), and using the expansion of the intensity distribution, in 
temperature units instead of intensity units, one can show g(x) — (x/2) coth(a;/2). 

The a22 quadrupole moment related to the cluster transverse motion (in a coordinate 
system in which z is along the line of sight) is 



(67) 



»22 




where <f> v is the orientation angle for vt on the plane of the sky. To obtain this result, 
note that in a coordinate system in which the z axis is aligned with the cluster's mo- 
tion, the quadrupole dependence of the radiation temperature is g(x)v 2 ([i 2 — 1/3) = 
g(x)v 2 (2/3)y/ '47r/5Y2,o, where (i is the cosine of the angle between the radiation direc- 
tion and the z direction. However, in a coordinate system in which the i axis is taken 
to be along the line of sight, (/x 2 - 1/3) = -(1/3)^4^/5^2,0 + y/2-K /\h{Y 2 , 2 + *2 -2)- 
Thus, the coefficient of Y 2>2 , the component of the quadrupole moment that gives rise to 
polarization in our direction, is only a fraction y/3/8 of the total quadrupole moment. 

Fig. illustrates the primordial and kinematic quadrupole contributions using all- 
sky maps of the expected polarization. In these plots, each polarization vector should 
be considered as a representation of the polarization towards a cluster at that loca- 
tion. The polarization pattern created by scattering of the primordial quadrupole is 
uniform and traces the underlying temperature quadrupole distribution. For the kine- 
matic quadrupole we assume, for illustrative purposes, a transverse velocity field with 
(fit) ~ 10~ 3 corresponding to a velocity of 300 km sec -1 . The polarization contribution 
due to the kinematic quadrupole, however, is random due to the fact that transverse ve- 
locities are uncorrelated( 2 ). This explains the randomization of the all-sky polarization 
map when a significant kinematic contribution is included. It should be noted, however, 
that the kinematic effect has been scaled by a factor of 100 to make it visible in Fig. [7| 



( 2 ) The correlation length of the velocity field (~ 60 Mpc) correlates velocities within regions 
of 1° when projected to a redshift of order unity. 
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Fig. 8. - Polarization power spectra due to the rescattering of the primordial and kinematic 
quadrupoles. We break total power in each of these cases into the 1- (lh; dotted lines) and 2- 
halo (2h; dashed lines) terms under the halo-based approach used here. While at large angular 
scales correlations between halos dominate, at small angular scales of order few arcminutes and 
below, contributions are dominated by the 1-halo term. For comparison, we also show primordial 
polarization power spectra for E and B-modes involving dominant scalar (E-mode) and tensor 
(B-mode) contributions respectively. The tensor contribution to B-modes due to inflationary 
gravitational waves (IGW) assumes an energy scale for inflation of 10 16 GeV. The long-dashed 
curve is the contribution to B-modes of polarization resulting from the cosmic shear conversion 
of power in E-modes, while the dot-dashed line labeled "Residual Lensing" represents the noise 
contribution after optimally subtracting the lensing contribution using higher-order statistics. 



Therefore, as shown in Fig. \7\ the primordial polarization dominates the total contri- 
bution even at high frequencies where the kinematic quadrupole is increased due to its 
spectral dependence. Also, the spectral dependence of the kinematic quadrupole contri- 
bution, g(x), gives a potential method to separate the two polarization effects |58|. This 
is similar to component separation suggestions in the literature as applied to temperature 
observations, such as the separation of the thermal SZ-effect from dominant primordial 
fluctuations 

5'4. CMB-induced Cluster Polarization and Dark Energy. - Secondary CMB polar- 
ization in the direction of distant galaxy cluster may provide an opportunity of study the 
evolution of growth in the late, dark energy dominated universe 05] , [5U], 100] ■ This will 
constrain the time-evolution of the dark energy equation of state and therefore might 
contribute to a better understanding of the properties and the nature of dark energy. 

The basic idea is to use galaxy clusters as tracers of the local temperature quadrupole 
and statistically detect its rms value. Imagine a future experiment measuring the polar- 
ization towards resolved clusters. For each cluster we also measure its redshift and the 
optical depth through the thermal SZ effect. We do this for a large sample of clusters, 
bin the resulting data into redshift intervals and average over all sky. If the bin size can 
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be chosen small enough this will allow us to reconstruct the rms temperature quadrupolc 
as a function of redshift. Multi- frequency observations will allow us to separate the pri- 
mordial quadrupole from the contaminant kinematic contribution with only a factor of 
~ 2 enhancement in the instrumental noise |49| . 

For the best-fit ACDM model with 70% dark energy, the quadrupole leads to a maxi- 
mum primordial polarization of Pp r i m ~ 4.9r /xK. Since the kinematic polarization scales 
as PKin = 0.27.g(a:)(/3 t /0.001) 2 T /xK, we expect the CMB-induced signal to be dominant 
(as illustrated in Fig. EJ. Factors that could make the primordial and the kinematic 
polarization more comparable are the frequency boost of the kinematic effect, a more 
optimistic estimate of (fit) and a low value of the CMB-quadrupole. However, even if 
the signals were of comparable magnitude the random orientations and characteristic 
frequency dependence of the kinematic polarization would allow a reliable extraction of 
the primordial CMB-quadrupole. 
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Fig. 9. - Top panel: Projected ISW quadrupole as a function of redshift. The solid error bars as- 
sume a reconstruction with clusters down to 10 14 Mq in an area of 10 4 deg 2 with an instrumental 
noise of 0.1 /xK. The dotted lines show the cosmic- variance for an all-sky reconstruction com- 
puted from the number of independent volumes sampled by clusters at each redshift bin |54|. 
Bottom left: Parameter errors from the projected ISW quadrupole measurements, assuming 
w — const. The small ellipse is for the case shown in top panel (with cosmic variance added 
in quadrature), while the two larger ellipses assume a factor of 3 and 10 increase in the instru- 
mental noise contribution, respectively. For comparison, we also show the constraints expected 
from SNAP. When the most optimistic polarization information is added, SNAP's constraints 
on w improve by a factor of 3. Bottom right: Same as bottom left, but for wo and w a of the 
parameterization 1x1(2) = Wo + w a z/(l + z) and assuming an additional prior <t(Qm) = 0.01. 

In we assessed the potential detectability of the ISW signature with cluster po- 
larization, while in 60 , the measurement of cosmological parameters related to the dark 
energy using cluster polarization data from, say, the planned CMBpol mission, was dis- 
cussed. The measurement of dark energy properties is aided by the fact that one probes 
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the redshift evolution of the ISW contribution through the growth rate of the gravi- 
tational potential. The latter provides the most sensitive probe of dark energy when 
compared to all other cosmological probes considered so far [HOj- This comes from the 
fact that the growth rate of the gravitational potential is directly proportional to the 
dark energy equation of state, while quantities such as the distance or the growth factor 
involve, at least, one integral of this quantity. 

Establishing the ISW effect and reconstructing its redshift evolution poses a signifi- 
cant experimental challenge. The advent of polarization sensitive bolometers, however, 
suggests that a reliable reconstruction of the ISW signature is within reach over the next 
decade. Given the strong dependence of the ISW effect on the background cosmology, 
cluster polarization can eventually be used as a powerful probe of the dark energy. A 
detailed analysis of the potential of this method for extracting information on the dark 
energy equation of state, w(z), is presented in [HO] and illustrated in Figure El 

5'5. Power Spectra of SZ Polarization. - In addition to statistical studies to recon- 
struct the primordial quadrupole from CMB polarization measurement in the direction 
of resolved galaxy clusters, in the limit that clusters are unresolved, power spectrum 
measurements would provide an approach to study the polarization signal. In general, 
a polarization map will consist of a measurement of the Stokes parameters Q(n) and 
U(h) as functions of position n on some patch of the sky. We can construct Fourier 
components Q(l) and U(V) by 

(68) X(l) = J d 2 h e lhil X(h) , 

where X = Q,U, and 1 is a vector in the plane of a region of the sky sufficiently small to 
be considered flat. Since the Stokes parameters Q and U depend on the choice of axes, 
we consider the rotationally invariant combinations, 

E(l) = cos(2<fr)Q(l) + sin(20i)t/(l) 

(69) B(l) = cos(20i)Z7(l) ~ sin(20i)Q(l) , 



where 4>\ is the angle between 1 and the chosen x-axis in the plane of the sky. We define 
the angular power spectra Cf E and Cf B from 

(70) (Y(l)Y(l')) = (2nfS D (l-l')CT Y , 



where Y = E,B, and the angle brackets denote an average over all realizations of the 
density field. 

If the radiation quadrupole moment is smooth over the region of sky we are consid- 
ering, then 

(E(\)E(l<)) = (B(l)B(l')} 
(71) =i((Q(l)Q(l')} + (£/(lMl')}) , 

while 



(72) 



(E(l)B(l')) = \ «Q(l)t/(I')) - (Q(1M1')» = 0, 
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where the latter equality is consistent with parity conservation. These results can be 
derived by noting that if the quadrupole moment is constant, then the orientation of the 
polarization is constant. If so, we may choose our axes on the sky so that U — 0. Then, 
E(\) oc cos(20i), and B(l) oc sin(2</>i), but when averaged over the orientation angle <f>\, 
we recover Eq. lf?T|) . 

We now proceed to calculate the power spectra induced by reionization. The polar- 
ization in direction n due to scattering from free electrons is an integral along the line of 
sight ED, 

(73) Q{n) - iU{n) = J d / T( g' v) a 22 ( v ) , 

where 77 is the comoving distance, (dT/dr/)(r]ii,ri) = OT"n e (r\xi, 77)0(77), 0(77) is the scale 
factor at a comoving distance 77, n e (776,77) is the free-electron density at direction n at 
distance 77, and <tt is the Thomson cross section. 

In Eq. I|73l) . 022(77) is the radiation quadrupole moment at distance 77. More precisely, 
022(77) is the coefficient of the spherical harmonic 122(0, </>) m a spherical-harmonic ex- 
pansion of the radiation intensity in a coordinate system in which the line of sight is the 
z direction. Note that we take 022(77) to be a function of distance only, and not direction, 
consistent with our assumption that the quadrupole is coherent over a large patch of 
the sky. Since we use Limber's approximation below, in which angular correlations are 
induced only by spatial separations at the same distance, the variation of 022(77) with 
distance can be included consistently. 

With the polarization written as a projection along the line of sight, Eq. 1)73(1. the 
angular power spectrum follows in the flat-sky limit from Limber's equation |21j . 

C EE — C BB 

where Ptt is the power spectrum of dr/drj, proportional to the power spectrum of the 
electron density n e , and the integral is taken up to the redshift z re ; of reionization using 
the comoving differential volume given by d 2 V/dfldz. 

In Fig- [HI w e summarize results related to the polarization power spectra from galaxy 
clusters. In the case of the primordial quadrupole, we replace |o22(^)| 2 in Eq. (|74() by 
its expectation value, Cf (z), the variance of the temperature quadrupole at redshift z 
(shown in Fig. EJ. In the case of the kinematic quadrupole, we replace |a22 (^) | 2 by the 
expectation value of equation (|67H , 

(75) <|a 22 | 2 ) = %9\*){vi) = ||.9 2 (*K 4 ms , 

where we have used (vf) = {(v% + v 2 ) 2 ) = (8/9> r 4 ms , since (v%) = (w 2 ) = w 2 ns /3, (v$) = 

( v t,) = 3( u2 ) 2 = u mis/3, an d (vl v l) = (v 2 )^ 2 .) = v f ms /9- We calculate the linear-theory 
rms peculiar velocity by integrating over the power spectrum, as in the case of the SZ 
kinetic temperature anisotropy calculation. According the halo-clustering model 
peculiar-velocity fields are correlated over large distances and the non-linear corrections 
to v^ ms are small. The resulting linear-theory rms kinematic quadrupole is also shown 
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in Fig. assuming g{x) = 1 as relevant for Rayleigh- Jeans (RJ) part of the frequency 
spectrum when x — > 0, and for v — 150 GHz and v — 220 GHz. 

In Fig. |H1 note that the secondary polarization discussed here contributes equally to 
E- and B-modes (cf. equation (|T4|> ). While the 1-halo term dominates at arcminutc 
angular scales and below, correlations between halos are important and determine the 
total effect due to secondary polarization at angular scales corresponding to a few degrees. 
The dependence of our results on the inclusion of the 2-halo term is consistent with the 
result obtained for the temperature power spectra from the kinetic SZ effect Jl], while 
it is inconsistent with the thermal SZ effect, where contributions are dominated by the 
I-halo term over the whole range of angular scales. The latter behavior is explained 
by the fact that the thermal SZ effect is highly dependent on the most massive halos, 
while the kinetic SZ effect, and the secondary polarization signals calculated here, are 
independent of the gas temperature and thus can depend on halos with a wider mass 
range. 

As shown in Fig. |H1 the secondary E-mode polarization is several orders of magnitude 
below the E polarization from the surface of last scattering. The secondary polarization is 
therefore unlikely to be a source of confusion when interpreting polarization contributions 
to E-modes. The amplitude of the primary effect in B-modes, due to gravitational waves, 
is highly uncertain and depends on the energy scale of inflation |fi2| . For illustration, we 
show in Fig. the inflationary gravitational wave (IGW) signal assuming an energy scale 
for inflation of E- ln n = 1 16 GeV, where the amplitude of the power spectrum scales as 
Ef nfi . At large angular scales the secondary polarization is several orders of magnitude 
below the peak of this hypothetical IGW polarization signal. If the energy scale of 
inflation is lowered considerably, say to E- m ^ < I0 15 GeV, then we might guess that the 
secondary polarization could ultimately constitute a background. 

As also shown in Fig. |S| however, there is a contribution to the B-mode power spec- 
trum that arises from conversion of the primary E-modes to B-modes by gravitational 
lensing |63|. and this is considerably larger than the secondary polarization from galaxy 
clusters. Moreover, we also show (the dot-dash curve) the contribution to the irre- 
ducible B-mode power spectrum that remains even after the lensing has been optimally 
subtracted with higher-order correlations ESI E23 EH • This residual lensing power 
spectrum is considerably larger than the polarization from reionization. Thus, the sec- 
ondary SZ polarization effects are unlikely to be a factor for either gravitational-lensing 
or gravitational-wave studies with B-modes. 

* * * 

AC thanks the organizers of the summer school, F. Mclchiorri and Y. Rcphacli, for 
an enjoyable and productive conference and for the invitation to present the research 
discussed in this review. We would like to thank Marc Kamionkowski and Dragan Huterer 
for important collaboration upon which parts of this review are based. 
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